Spatial state Stokes-operator squeezing and entanglement for optical beams 
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The transverse spatial attributes of an optical beam can be decomposed into the position, mo- 
mentum and orbital angular momentum observables. The position and momentum of a beam is 
directly related to the quadrature amplitudes, whilst the orbital angular momentum is related to 
the polarization and spin variables. In this paper, we study the quantum properties of these spa- 
tial variables, using a representation in the Stokes-operator basis. We propose a spatial detection 
scheme to measure all three spatial variables and consequently, propose a scheme for the generation 
of spatial Stokes operator squeezing and entanglement. 

PACS numbers: 42.50, 42.30 



I. INTRODUCTION 

Squeezed and entangled bright optical fields are es- 
sential resources in the continuous variable quantum op- 
tics and quantum information communities To date, 
the vast majority of research has been focused on fields 
which exhibit non-classical features on their amplitude 
and phase quadratures. Recently, however, a number of 
papers have been published on squeezing and entangle- 
ment of other field variables, and in particular the polar- 
ization [U and spatial structure [!, 0, H, II] • These vari- 
ables can be directly related to momentum components of 
the field. The quadrature amplitudes are associated with 
the transverse linear momentum of the field, whilst the 
polarization states and spatial states considered to date 
are respectively associated to the spin angular momen- 
tum and the transverse angular momentum. We immedi- 
ately see why interaction of a polarization squeezed field 
with an atomic ensemble can yield atomic spin squeezing 
, and why spatially entangled light can be used to test 
the EPR paradox with position and momentum variables 
[f| as was originally proposed by Einstein, Podolsky and 
Rosen |[. 

The correspondence of non-classical polarization and 
spatial states with non-classical momentum states imme- 
diately raises the question of whether non-classical orbital 
angular momentum states can also be generated. Such 
states have been under investigation for some time in the 
discrete variable regime. Techniques have been estab- 
lished to detect the orbital angular momentum proper- 
ties of single photons [1, [l(| EH; and discrete- variable 
multi-dimensional entanglement between orbital angular 
momentum states has been proposed [l2l [l3j and demon- 
strated However, to date, non-classical orbital angu- 
lar momentum states have not been investigated for the 
continuous variable regime, that is relevant in many ap- 
plications [H [M EE UM, El IM EH . We investigate these 
continuous variable orbital angular momentum quantum 
states, proposing techniques for both generation and de- 
tection. Many applications have been proposed for op- 
tical orbital angular momentum in the classical domain, 



ranging from the generation of counter-rotating superpo- 
sitions in Bose- Einstein condensates [12, HH , and transfer 
of orbital angular momentum from an atomic ensemble 
to a light field [24] , to optical angular momentum trans- 
fer to trapped particles [Ifl [2(| [27| • It should be noted 
that light with orbital angular momentum has also been 
applied to achieve super-resolution imaging of molecules 
and proteins in biological systems, via techniques such as 
stimulated emission depletion (STED) 28]. Non-classical 
orbital angular momentum states offer the prospect to 
improve these processes, and to generate non-classical or- 
bital angular momentum states in macroscopic physical 
systems. 

Ref. [5[ proposed EPR entanglement of the position- 
momentum observables of optical fields, and this has re- 
cently been experimentally observed [6|. In these ref- 
erences, entanglement was present only for one trans- 
verse axis of the beam, and analysis only required con- 
sideration of the TEMoo and TEM r modes. One of the 
strengths of entanglement in the spatial domain is access 
to higher dimensional spaces. Here we extend the work 
of Refs. [1,0 to include both transverse beam axes, in a 
formalism easily extended to higher order TEM modes. 
In the process, orbital angular momentum is introduced 
as an entanglement variable in addition to beam position 
and momentum. These quantum variables can be repre- 
sented on a spatial version of the Poincare sphere com- 
monly used for analysis of polarization quantum states. 
We propose a scheme to measure their signal and noise 
properties, via a spatial Stokes detection scheme; and 
finally introduce a scheme to generate spatial Stokes- 
operator squeezing and entanglement. 



II. SPATIAL STOKES OPERATORS 

Spatial quantum states exist in an infinite dimensional 
Hilbert space, which may be conveniently expanded in a 
basis of TEM P9 modes. After such an expansion the pos- 
itive frequency part of the electric field operator, £ + (r), 
can be written explicitly in terms of the double-subscript 
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where a p<3 and u pq (r) are the photon annihilation opera- 
tor and the normalized transverse beam amplitude func- 
tion, respectively, associated with the TEM pg mode. The 
photon annihilation operator be formally defined 

with the projection of the transverse beam amplitude 
function with the positive frequency part of the electric 
field operator, given by 



dxdy£ + (x,y)u pq (x,y). 



(2) 



Within this basis, it is natural to consider pairs of modes 
with amplitude function that correspond to a physical 
7r/2 rotation in the transverse plane. That is, modes of 
the form TEM pg and TEM 9P , with p and q interchanged 
(p Pairs of this form are naturally described by a 

set of spatial Stokes operators S'f p ' 9 ^ where i g {0, 1, 2, 3}, 
that can be represented by a Poincare sphere |29l| in direct 
analogy to polarization Stokes operators [30|, as shown 
in Fig. ED 




FIG. 1: (a) Poincare sphere representation based on the spa- 
tial modes TEMio and TEMoi . Si is the Stokes variable for 
the uoi(r) and Mio(f) modes, S2 is the Stokes variable for 
the «oi ( r ) and wfo (r) modes, and S3 is the Stokes vari- 
able for the Wgl(r) and Uoi(r) modes, (b) Poincare sphere 
representation for the spatial Stokes operators. The shaded 
area indicates the quantum noise associated with the mean 
amplitude of the Stokes operator. 



where 9 is the phase difference between the TEM p? and 
TEM 9P modes. As is the case with polarization, the spa- 
tial Stokes operators obey commutation relations. By 
using the commutation relations of the photon annihi- 
lation and creation operators, [d mn ,ajk] — S m jS n k, the 
Stokes-operator commutation relations can be found to 
be 



Si^ Sj - — - 2z6 



i j k $k 



(4) 



where e^-fc is the Levi-Civita tensor. 

Each Stokes operator characterizes a different spatial 
property of the field. Sq represents the beam inten- 
sity, whilst the Stokes vector (Si, S2, S3) characterizes its 
transverse momentum. Si and £2 respectively quantify 
the relative proportions of horizontal to vertical trans- 
verse momentum, and 45° to —45° diagonal transverse 
momentum; whilst S3 weights left to right orbital angu- 
lar momentum. 

The spatial Poincare sphere can be fully spanned by 
spatially overlapping two orthogonal TEM pg and TEM gp 
modes. The diagonal modes, for example, are given by 



(r)- 



(r) and?C (r) = u qp (r) + u pq (r) , 



a pq ( r ) — wqpV I ~ "p<j 

and respectively yield Stokes vectors oriented along the 
negative and positive S% axis of the Poincare sphere. The 
Laguerre-Gauss modes LGoq with orbital angular mo- 
mentum q and —q, are given by u^ q (r) — Uo q (r)+iu q o(r) 
andug '(r) = ua q (r)—iu q0 (r), yielding Stokes vectors ori- 
ented along the S3 axis. For example, the intensity dis- 
tributions for the diagonal and Laguerre-Gauss modes, 
generated from the combinations of TEM10 and TEM01 
modes, are shown in Fig. [2j 
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Using the definition of the classical Stokes parameters 
[2{|, corresponding quantum mechanical spatial Stokes 
operators are defined as 
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FIG. 2: Intensity distribution for the modes given by (A) 
"of (r) = uoi(r) +uio(r), (B) u\l (r) = Moi(r) -uio(r), (C) 
u oi( r ) = u 01 (r)+iu 10 (r), and (D) u^(r) = uoi(r)-iuw(r). 



Since the diagonal and orbital modes can be generated 
purely from two orthogonal TEM pg and TEM gp modes, 



only by changing the phase between the modes, these sets 
of modes obey the SU(2) group properties. At this point, 
we restrict our analysis to the TEMi and TEM i modes, 
for illustrative purposes. In principle, our analysis is valid 
for all orthogonal TEM P9 and TEM gp modes. 



III. SPATIAL STOKES DETECTION 



Now that the spatial Stokes operators have been de- 
fined, we consider a system for their detection. Analo- 
gous to polarization Stokes operator detection, the spa- 
tial Stokes detection requires two separate photodiodes, 
a modal phase shifter, and a mode separator. 

A modal phase shifter is a device that introduces a 
relative phase shift between the two orthogonal TEMio 
and TEMoi modes. The modal phase shifter can be con- 
structed using a pair of cylindrical lenses with variable 
lens separation, as shown in Fig. [3J (f). The pairing of 
cylindrical lenses introduces an astigmatic Gouy phase 
shift between the TEMi and TEM i modes. In order to 
introduce 7r and 7r/2 modal phase shift, the lens separa- 
tion is given by 2/ and y/2f, respectively, where / is the 
focal length of the cylindrical lens. Ref. [2l[ contains a 
detailed analysis of the tt and 7r/2 modal phase shifters. 

A mode separator (MS) is a device that separates the 
two orthogonal TEMio and TEMoi modes and can be 
constructed using an asymmetric Mach-Zehnder inter- 
ferometer |':jrO • as shown in Fig. [3J (e). Due to the odd 
and even numbers of reflections in each interferometer 
arm, different interference conditions for the TEMoi and 
TEMio modes are present. The resulting outputs from 
the interferometer is a separation of even and odd (de- 
fined along one spatial axis) spatial modes [3ll |. 

Measurements of the total signal and noise, which cor- 
respond to a measurement of So , are given by the sum of 
the photocurrents as shown in Fig.[3j(a). Measurement of 
Si involves taking the subtraction between the photocur- 
rent outputs corresponding to mode components TEMio 
and TEMoi, as shown in Fig.[3J(b). Measurement of 52 
involves subtraction of the diagonal mode components 
and is obtained by phase shifting one mode by tt with 
respect to the other and then taking the subtraction of 
the photocurrent signals, given in Fig. [3J (c). S3 is mea- 
sured by decomposing the Laguerre-Gauss mode into its 
TEMio and TEMoi modes using tt and ir/2 modal phase 
shifters, as shown in Fig. [3J (d). 

The photon annihilation operators in Eq. (J3J can be 



written in the form a pq = a pq - 



-8a pq , where a pq describes 



the mean amplitude part and Sa pq is the quantum noise 
operator. Using the linearized formalism, where second 
order terms in the fluctuation operator are neglected (i.e. 
\ a pq\ 2 ^ K^m)I)j tne mean amplitudes of the Stokes 
operators in Eq. ©, in terms of the TEMio and TEM i 
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FIG. 3: Measurements of (a) So, (b) Si, (c) S 2 and (d) S3. 
(e) An example of a mode separator (MS) is an asymmetric 
Mach-Zehnder interferometer, (f) The n and tt/2 modal phase 
shifters could be constructed using a pair of cylindrical lenses 
with lens separation given by 2/ and \/2f, respectively. / is 
the focal length of the cylindrical lenses, M is a mirror, 50:50 
BS is a symmetric non-polarizing beam-spliter. 



modes, are therefore given by 
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(Si) 
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a 1Q - a 
2aioQ!oi cos 6 
2aioQfoi sin 9 



(5) 



where aoi and aw are the mean amplitude terms cor- 
responding to modes Moi( r ) and u\o(r), respectively. 
(So) = N is the total mean number of photons, (S%) is the 
difference in the mean number of photons in the uio^) 
and U01 (r) modes, (S2) is the difference in the mean num- 
ber of photons in the Uqi° ( r ) an( i u to° ( r ) modes and (S3) 
is the difference in the mean number of photons in the 
v,Qi(r) and Uoi ( r ) m °des. 

The Stokes operators of Eq. (J3J can be expanded in 
terms of quadrature operators with the general form 



I„, and their variances are then 
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given in general by 

((SS ) 2 ) = a 2 o<(«i ) 2 )+<((« a + 01 ) 2 > 

+2a 10 a i(^i «a + 01 > 
((^) 2 ) = a 2 1Q ((SX+J 2 )+a 2 01 ((6X+ Qi ) 2 ) 

-2a w a Q1 (SX+ o SX+J 
((SS 2 ) 2 ) = ai ((5Xjj 2 ) + aU(6X e a J*) 

+2a w a 01 (5X- o e i 6X e a J 
((6S 3 ) 2 ) = al ((6X^) 2 ) +aU(8X e a -%f) 

+2a 10 a 01 (SXa^ 2 SXa 10 2 ) 

(6) 

where = Xf =0 and XT' = if" 1 '' 2 are respec- 

tively the amplitude and phase quadrature operators of 
the TEM M mode. 

Non-classical optical orbital angular momentum states 
can be constructed by spatially overlapping a set of or- 
thogonal non-classical TEM p(3 fields, in analogy to the 
work of Refs. 0, for transverse spatial entanglement. 
In such a scenario, no correlations should exist between 
the quadratures of the different input fields, such that 
(SX^JX^) = 0, V-j>i,c/>2}. Making this assumption, 
which we will adopt henceforth, the variances of the 
Stokes operators are simplified to 

((SS ) 2 ) = aUiSX+^+aUiSX+J 2 ) 
((6S,) 2 ) = a 2 10 ((SX+J 2 )+a 2 Q1 ((6X+ m ) 2 ) 
((SS 2 ) 2 ) = aUiSX^+aUiSX^J 2 ) 
((SS 3 ) 2 ) = a 2 ((SX-^) 2 ) +a 2 )1 ((5X e a -t) 2 ). 

(7) 

IV. SPATIAL STOKES SQUEEZING 

An examination of Eq. ([7]) shows that simultaneous 
squeezing of at least two spatial Stokes operators is pos- 
sible, using two quadrature squeezed input fields. Hence, 
a spatial Stokes squeezed state can be used to enhance 
relative measurements of momentum variables along mul- 
tiple axes. This is best illustrated by considering a few 
examples. 

Simultaneous squeezing of 5*0, Si and S 2 can be 
achieved through the in-phase (6 = 0) spatial over- 
lap of amplitude squeezed TEMio and TEMoi modes 
(i.e. {{5X+ Q f) < 1 and ((SX+J*) < 1). This spatial 
squeezed state could be combined with a bright TEM o 
beam to allow enhanced measurements of the transverse 
momentum of the TEMoo beam along any axis on the 
Poincarc sphere within the plane formed by Si and 52. 
Simultaneous squeezing is exhibited along the horizon- 
tal/vertical axis and the diagonal/anti-diagonal axis. 

Squeezing of So, Si and S3 can be achieved by overlap- 
ping amplitude squeezed TEM10 and TEMqi modes with 



phase difference 9 — ir/2. The combination of this spatial 
squeezed beam with a bright TEMoo beam would simul- 
taneously enable sub-shot noise relative measurements of 
transverse momentum between the horizontal and verti- 
cal well as left- and right-handed orbital angular 
momentum of a TEMoo beam 

To achieve quantum enhanced measurements of trans- 
verse momentum between the diagonal and anti-diagonal 
axes, as well as between left- and right-handed orbital 
angular momentum of a TEMoo beam, would require 
squeezing of S 2 and S3. This can be achieved by over- 
lapping quadrature squeezed TEM10 and TEM01 modes, 
at squeezed quadrature angle of 7r/4 (i.e. ((<5A a 4 10 ) 2 ) < 1 
and ((5X? 01 ) 2 ) < 1) and 6 = tt/4. 

Phase squeezed TEM10 and TEM01 modes (i.e. 
((Si-J 2 > < I and {{6X-J 2 ) < I) will yield squeez- 
ing of ((SS 2 ) 2 ) only. 

V. SPATIAL STOKES ENTANGLEMENT 

Spatial Stokes entanglement can be generated as is 
shown in Fig. [U Two squeezed TEM10 beams, labeled 
respectively by the subscripts x and y, are combined with 
a 7r/2 relative phase shift on a 50:50 beam-splitter. The 
beams after the beam-splitter, also in TEM 10 modes, ex- 
hibit the usual quadrature entanglement. Two spatial 
mode combiners, each consisting of a spatial mode sep- 
arator as shown in Fig. [3] (e) in reverse, are then used 
to combine the quadrature entangled beams with bright 
coherent TEM01 beams, labelled here <Zoi !X and aoi, a , re- 
spectively. As we will show here, the resulting output 
beams are entangled in the spatial Stokes operator ba- 
sis. Experimental verification of the spatial Stokes en- 
tanglement can be performed using the detection scheme 
described in the preceding sections. 




FIG. 4: Scheme to generate and characterize spatial Stokes 
entanglement. 50:50 BS is a symmetric non-polarizing beam- 
splitter, 7r and 7r/2 are modal phase shifters. 

To verify the existence of spatial Stokes entanglement 
we will use the generalized version of the Duan insepa- 
rability criterion [32j given in Ref. [HI]. This criterion 
provides a sufficient condition for entanglement and is 
given by 

A 2 x±y A + A 2 x±y B<2\[SA,5B}\ (8) 



5 



where A and B are two general observables, and 
A 2 x±y 6 = mm((SO x ± (SOy) 2 ). The degree of insepa- 
rability 1(A, B) is then given by [33j 



A x± a ^ + A 2 x±y B 

2\[SA,SB}\ ' 



(9) 



where X(.A, -B) < 1 indicates an inseparable state. 

Assuming for simplicity and symmetry (i.e. both 
beams x and y are identically interchangeable) that 
uoi.x = «oi,y = «oi and ai ,x = aio,y = «io, and us- 
ing Eqs. f4} and ([5]), we find the inseparability criteria 
between spatial Stokes operators given by 



X(S h S 2 ) = 
I(S 3 ,Sx) = 
I(S 2 ,S 3 ) = 



A x ± y Sl + A 2 ±y S 2 

8|aioaoi sin0| 

A x ± y S\ + &l± y S3 

8|aioaoi cos $| 

<^l± y S 2 + ^l±y S 3 

4 Ko- a oil 



(10) 



The correspondence between spatial Stokes and 
quadrature entanglement becomes obvious when we ex- 
press the spatial Stokes operator conditional variance 
A 2 ±y Si in terms of quadrature operators. Making the 
assumption that aio <IC aoi and using Eq. ([5]), gives 

= a 2 al Al ±y X+ oi 

^■x±y^2 = a 01^x±y^a 10 

A 2 c _ n 1 A 2 V d ~i 
^■x±y C3 3 — a 01 lA x±y jL a 10 ■ 

Substituting these expressions into Eq. (fTU|) . gives 



(11) 



T(S U S 2 ) 



A 2 v-+ , a 2 y6 

8ai |sin6»| V *±« a <» ^ x± v a ™ 



«01 



8ai O |cos0| V ^ 



/ A 2 y+ , a 2 v- 6 "! 



(12) 



I(5 , 2,S , 3) — ^{A x ±yX e aw + A 2 x±y X aia 



(13) 



(14) 



An inspection of Eqs. (fT2f to (fT4|) show that in the 
limit considered here only spatial Stokes entanglement 
between S 2 and S3 can be realistically attained. This is 
because the inseparability criteria of Eqs. (|12[) and (fT5|) 
scale with respect to aoi/aio- Since aoi/aio 3> 1, the 
correlation in X+ Qi and X® w has to be significantly re- 
duced below one, in order to maintain inseparability. On 
the other hand, we see that after setting 9 = 0, Eq. (fl4|) 
reduces to 



1{S 2 , S 3 ) - - (A£± y X£ 10 + A x±y X au 



(15) 



So that quadrature entanglement between the TEM10 
modes is transformed directly into spatial Stokes entan- 
glement between the S 2 and S3 spatial Stokes operators. 



VI. CONCLUSION 



We have identified the relevant spatial modes for op- 
tical beam position, momentum and orbital angular mo- 
mentum. We formalized the quantum properties of these 
observables using the Stokes-operator formalism and pre- 
sented a graphical representation of these variables using 
the Poincare sphere. A spatial Stokes detection scheme 
was described and schemes to generate spatial Stokes op- 
erator squeezing and entanglement were proposed. 
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